Gauge invariance in two-particle scattering^ 



B. Borasoy2«, P. C. Bmns^", U.-G. MeiBner^"'^ R. NiBler^' 

"Helmholtz-Institut fiir Strahlen- und Kernphysik (Theorie) 
Universitiit Bonn, NuBallee 14-16, D-53115 Bonn, Germany 

''Institut fiir Kernphysik (Theorie), Forschungszentrum Jiihch 
D-52425 Jiihch, Germany 



Abstract 

It is shown how gauge invariance is obtained for the coupling of a photon to a two-body 
state described by the solution of the Bethe-Salpeter equation. This is illustrated both 
for a complex scalar field theory and for interaction kernels derived from chiral effective 
Lagrangians. 
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1 Introduction 



Chiral perturbation theory provides an appropriate framework for studying hadronic processes 
at low energies [1] . In chiral perturbation theory (ChPT) the most general effective Lagrangian 
incorporating the same symmetries and symmetry breaking patterns as the underlying theory, 
QCD, is formulated in terms of the relevant degrees of freedom, i.e. mesons and baryons. 
Moreover, external fields representing, e.g., external photons are included in a gauge invariant 
fashion. The Green's functions are then expanded perturbatively in powers of Goldstone boson 
masses and small three-momenta. By employing a regularization scheme which respects chiral 
symmetry, gauge invariance is maintained at every order in the loop expansion of ChPT. 

However, the systematic loop expansion involves a characteristic scale ^ 47rF7r ~ 1.2 GeV 
at which the chiral series is expected to break down and the limitation to very low-energy 
processes is even enhanced in the vicinity of resonances. The appearance of resonances in 
certain channels constitutes a major problem for the loop expansion of ChPT, since a resonance 
cannot be reproduced at any given order of the chiral series. Nevertheless, at low energies the 
contribution from such resonances is encoded in the numerical values of certain low-energy 
constants (frequently called resonance saturation). 

Recently, considerable effort has been undertaken to combine the effective chiral Lagrangian 
approach with non-perturbative methods, both in the meson-baryon sector [2-4] and in the 
purely mesonic sector [5]. The combination with non-perturbative schemes has made it possible 
to go to energies beyond and to generate resonances dynamically (giving up, however, 
certain aspects of the rigorous framework constituted by ChPT). Two prominent examples in 
the baryonic sector are the A(1405) and the 5'ii(1535). The first one is an s-wave resonance just 
below the K~p threshold and dominates the interaction of the KN system, while the 5'ii(1535) 
was identified in [2] as a quasi-bound KK-KT^ state. 

Such chiral unitary approaches have been extended to photo- and clectroproduction pro- 
cesses of mesons on baryons, see e.g. [6-8]. In these coupled channel models the initial photon 
scatters with the incoming baryon into a meson-baryon pair which in turn rescatters (elastically 
or inelastically) an arbitrary number of times. The two-body final state interactions are taken 
into account in a coupled-channels Bethe-Salpeter equation (BSE) or - in the non-rclativistic 
framework - Lippmann-Schwinger equation. The coupling of the incoming photon to other pos- 
sible vertices is omitted, and although these approaches appear to describe the available data 
well, the issue of gauge invariance is not discussed in these chiral unitary approaches. On the 
other hand, a method to obtain conservation of the electromagnetic current of a two-nucleon 
system is presented in [9] and extended to a resonance model for pion photoproduction in [10]. 
Alternatively, the so-called "gauging of equations" method has been developed in [11,12] to 
incorporate an external electromagnetic field in the integral equation of a few-body system 
in a gauge invariant fashion. Further investigations of gauge invariance in pion photopro- 
duction within t:N models can be found in [13, 14]. Gauge invariance is also of interest in 
coupled-channels approaches in the mesonic sector, e.g in radiative 0, p decays [15,16] and in 
the anomalous decays 77,77' 77,7r"'"7r~7 [17,18]. For related recent work, see also [19]. 

The purpose of the present work is to illustrate how gauge invariance can be maintained 
when an external photon couples to a two-particle state described by the BSE within the chiral 
effective framework. We will start in the next section by first outlining the procedure for a 
simple scalar field theory. With the insights gained from this example we can then address 
gauge invariance in meson-baryon scattering processes with chiral effective Lagrangians. In 
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Sec. ISlwe discuss the case with the interaction kernel of the BSE derived from the leading order 
Lagrangian - the Weinberg- Tomozawa term. The extension to driving terms from higher chiral 
orders is presented in Sec. |3] In Section El it is shown that the corresponding amplitudes also 
satisfy unitarity constraints. Our conclusions and outlook are presented in Sec. El 



2 Complex scalar fields 

In this section, we will demonstrate in a simple field theory with complex scalar fields that the 
coupling of an external photon to a two-body state, which corresponds to the solution of the 
BSE, is gauge invariant. To this end, consider the (normal-ordered) Lagrangian for complex 
fields (j) and ip 

C = 9^0*9^0 - m^(f)*(j) + 8^-^*8"'^ - M^iIj*iIj - g{(f)*(j)) (ifj*-^) , (1) 

with masses m and M, respectively. For small values of the coupling constant g, the scattering 
process (p{pi)ip{p2) 0(P3)V'(P4) be calculated perturbatively. For general values of g, 
however, and if one is interested in bound states, one must resort to non-perturbative techniques 
such as the BSE. The Bethe-Salpeter equation for the scattering matrix T of the two-particle 
scattering process can be written as 

Tis) = g + gGis)T{s) , (2) 

where s = = {pi + P2Y = {ps + P^)"^ and G is the scalar loop integral 

G{p^)=tJ^A4l)A4p + l) (3) 



utilizing the short-hand notation 



i2n) 
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(4) 



and the propagators 



The solution of the BSE is given by 

^(") = l-gG{s) = ^ + ^^^'^^ + gGis)gGis)g + ... (6) 

which amounts to the summation of an infinite series of s-channel loops, the so-called bubble 
chain (or bubble sum) . From inversion of Eq. © it immediately follows that 

ImT^^ = -ImG = ^^^(s-(m + M)2) , (7) 

where the second relation is deduced from unitarity of the amplitude T and q^^ is the three- 
momentum in the center-of-mass frame. The explicit calculation of the scalar loop integral G, 
e.g. in dimensional regularization, indeed confirms the above relation. 
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Let us now turn to the coupling of an external photon field to the solution of the BSE. By 
minimal substitution 

dn ^ 1^ ^ + ie^Afj, for , 

d^i^Vi^^ + ie^A^ for V' , (8) 

where is the photon field and (e^) the charge of the meson (p H'), one obtains a locally 
gauge invariant Lagrangian. The coupling of the photon with incoming momentum k to the 
four external legs of a bubble chain leads to the amphtudes 

T^'^e^T{s')A{p^ + k) {2p^ + kY , 

= T{s') A{p2 + k) {2p2 + kf , 
T^ = e^{2p^-kYA{p^-k)T{s) , 

= ev,(2p4 - kY A(P4 - k) T{s) , (9) 

where s' = {p + k)"^. Multiplying these contributions with the four-momentum of the photon, 
k/^, and setting the external legs on-shell yields 

(e^ + e^) [T{s') - T{s)] , (10) 

which in general does not vanish. This underlines that, in order to achieve gauge invariance, it is 
not sufficient to couple the photon only to external legs. One rather has to include contributions 
which arise due to the coupling of the photon to intermediate states within the bubble chain, 
leading to the additional contributions 

T,^ = te^T{s') ^ A^{1 + k) {21 + kf A^{1) A^{p - I) T{s) , 

= te^T{s') ^ A^l + k) {21 + kY A^{1) A^{p - I) T{s) . (11) 

By employing the Ward-Takahashi identities 

k,{2l + kY = A-}^{1 + k)- A-}^{1) (12) 
it is straightforward to show that 

k, {T^ + Ti:) = (e^ + e^) T{s') [G{s) - G{s')] T{s) . (13) 
The last expression can be rewritten by making use of the BSE 

T{s')[G{s) - G{s')]T{s) = T{s') [g-'T{s) - 1] - [T{s')g-' - 1] T{s) = T{s) - T{s'). (14) 
Adding up all contributions we arrive at 

6 

1=1 

which confirms the gauge invariance of the photon coupling to the bubble chain. 
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At the same time, insertion of the photon coupUng at all possible places in the bubble chain 
guarantees unitarity of the scattering matrix up to radiative corrections of order 0{e^). To this 
end, we remark that in the transition (pip'-f — > (pip we can restrict ourselves to the states \(f)ip) 
and {(pt/j'-f) . Unitarity of the iS-matrix, S = 1 — iT, implies then 

(</)^|r-rt|</,^7) = -tJ |(0^|rt|</,y)(</,y|r|0^7) + (0V|t^I</>W)(</>W|T|</)^7)} , 

(16) 

where Jp^ denotes the phase space integral for the set of intermediate states and I^V't')- 

We have introduced a superscript for the intermediate particles with running three-momenta, 
(j)', ip' and 7', in order to distinguish them from the external particles 0, ip and 7 with fixed 
momenta. Note that the last matrix element {(j)'ip''y'\T\(j)ip'j) contains a disconnected piece of 
order 0{e^) in which the photon does not couple to the bubble chain and thus appears in the 
0{e) part of the unitarity relation. The remaining connected diagrams of this matrix element 
which contain the coupling of the photon to the bubble chain are of order 0{e^) and will be 
neglected in the following. By making use of the symmetry {i\T\j) = {j\T\i) due to time 
reversal invariance, Eq. |T6|l can be rewritten as 

-2im(0^|r|0^7)= / {(0^|T|0y)*(0y|r|0^7) + (0^|T|0y7)*(0y7|r|0v^7)| , (ir) 

Jps ^ ^ 

where now the phase space integral applies only to the particles (p' and ip'. The last equation 
represents the Cutkosky cutting rules [20]. At the diagrammatic level this amounts to cutting 
a pair of (p and ip propagators at all possible places in the bubble chain (keeping in mind that 
the photon is merely an external particle). The two terms on the right hand side represent the 
two possibilities for the photon to couple to the bubble chain before or after the cut. Since 
these are the only possible cuts in the bubble chain leading to imaginary pieces in the relevant 
kinematic region, one verifies Eq. ()17|) and hence unitarity of the 5-matrix up to radiative 
corrections. However, if the photon couples to a propagator, there is also the possibility to cut 
in the corresponding diagram both propagators which are directly connected to the photon. 
For a photon with fc^ < 4min(m^,M^) which is the case both for physical photons and for 
virtual photons from electron scattering, these cuts do not yield imaginary values and can be 
safely omitted here. 

Having convinced ourselves that it is possible to obtain gauge invariant and (up to radiative 
corrections) unitary amplitudes by taking into account the coupling of the photon to scalar 
fields in all possible ways in the bubble chain, we can now continue by applying this procedure 
to the slightly more complicated case of the chiral effective meson-baryon Lagrangian. 

3 Weinberg- Tomozawa term 

The chiral effective Lagrangian describing the interactions between the octet of Goldstone 
bosons (vr, K, 77) and the ground state baryon octet (A^, A, E, S) is given at leading chiral order 
by 

4b = ^{Bl,[D'', B]) - mo{BB) + . . . , (18) 

where (. . .) denotes the trace in flavor space. The 3x3 matrix B collects the ground state 
baryon octet and mo is the common baryon octet mass in the chiral limit. We have only 
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displayed the terms relevant for the present investigation and omitted two operators which 
contain the axial vector couplings of the mesons to the baryons. In the present investigation, 
we restrict ourselves to interaction kernels of the BSE given by contact interactions. The axial 
vector couplings of the mesons to the baryons could in principle contribute via direct and 
crossed Born terms, but the crossed Born term corresponds to three-body intermediate states 
which are beyond the scope of this work and we neglect the Born terms throughout. (The 
inclusion of Born terms in the interaction kernel is deferred to future work [21].) In fact, many 
coupled-channels approaches only take into account the contact interaction originating from 
the Lagrangian in ()18|). see e.g. [22] and references therein. 
The covariant derivative of the baryon field is given by 

[D,,B]=d,B + [V^,B] (19) 

with the chiral connection 

= ^ - ^ [y^^^i^u + uv^iS^ , (20) 

and = — e^^Q, where Q = |diag(2, —1, —1) is the quark charge matrix. The pseudoscalar 
meson octet (j) is arranged in a matrix valued field 

[/(</>) = u'i<f^) = exp (^V2z|^ , (21) 

with F the pseudoscalar decay constant in the chiral limit. Expansion of the chiral connection in 
the meson fields yields at leading order a (p'^BB contact interaction, the Weinberg- Tomozawa 
term, which we choose to be the driving term for the BSE in this section. 
The mesonic piece of the Lagrangian at leading chiral order is given by [1] 

4'^ = ^(V,f/W) + ^(x+), (22) 

where x+ = 2BQ{u'^Aiu^ + uM.u) describes explicit chiral symmetry breaking via the quark 
mass matrix M. = diag (m^, m^, m^), and Bq = ~ {0 \ qq |0) /F^ represents the order parameter 
of spontaneously broken chiral symmetry. The covariant derivative of the meson fields is given 
by (neglecting external axial-vector fields) 

V^t/ = d^U - iv^U + lUv^ . (23) 

In the Bethe-Salpeter formalism we choose to work with the propagators 



Sa{p) = , (24) 

with flavor indices i, a and physical meson and baryon masses. Mi and rria, respectively. How- 
ever, the following calculations are valid for all propagators satisfying the Ward-Takahashi 
identities 

k''V;l;{p + k,k) = A-\p + k)- A-\p) , 

k'^V^ip + k,k) = S-\p + k)~ S-\p) , (25) 
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where V^f (V^) are the corresponding 70^ {^BB) three-point functions. 

In the presence of a general interaction kernel A and coupled channels consisting of a meson- 
baryon pair the BSE for the process (f){qi)B{pi) —>■ (f){qf)B{pf) generalizes to 

I 

= Afi{p; qf, qi)+iY^ / Afi{p; qf, k) Aj{k) Saip + k) Tu^p; k, qi) (26) 
I -^^ 

with p = Pi + qi = Pf + qj and / = {(f)j,Ba} the channels which couple both to the initial and 
final state, i and /. Note that we have replaced the common mass of the ground state baryon 
octet, mo, by the physical baryon masses m^. This is consistent with the chiral order of the 
interaction kernel derived from the Weinberg- Tomozawa term and, in particular, produces the 
unitarity cuts at the physical thresholds. 




ABC 

\ / 



D E 

Figure 1: Tree diagrams for the process 705 —>■ (pB. Solid, dashed and wavy lines correspond 
to baryons, mesons and photons, respectively. The square denotes the vertex from the leading 
order Lagrangian. 

After setting up the formalism we first calculate the tree level contributions for the coupling 
of a photon to meson-baryon scattering. The pertinent Feynman diagrams for the process 
j{k)(j)i{qi)Ba{pi) 4>j{qf)Bb{pf) are depicted in Figure [H In addition to the coupling of the 
photon to the propagators the chiral connection in Eq. (fTHj) gives rise to a '-ycjP'BB vertex. 
Fig. EE. 
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The tree contributions to the transition amphtude read 

T^'^^'^''^' = "il^ { + if) ^'^(p^ + ^) ^^Q" + i^Q's^Pf - k) ii, + if) 

Hi, + + Mli + k) [2qi + k]^Q' 
Hi, + if-^) ^Mf - k) [2g/ - k],Q^ 

-1,{Q'+Q')\ {X^X'WXX']) , (27) 



where Q"A" = [Q, A"] (no summation over a) is the charge of the particle a in units of e and 
the A* are the generators of the SU(3) Lie- Algebra in the physical basis. By multiplying the 
tree contributions in Eq. (P7j) with gauge invariance is easily verified, if the momenta of the 
particles are put on-shell. 

In order to prove gauge invariance for the coupling of the photon to the bubble chain, it is 
convenient to consider first the diagrams presented in Fig. |21 with the pertinent contributions 
given by (a, 6, c, d denote baryon flavor indices, whereas z, j, m, n represent meson flavors): 

Fig. |2tA: 

4''"' = ^ 1,{Q' + Q') (A^niA^^ A'], A^> , (28) 

Fig. EB: 

// " <'^™5,(p - /) A^(/) T^-'"(p; /, q,) , (29) 

Fig.ET: 



Fig. EP: 



I j Af^^S^ip - /) A,^(/) r^™'»(p; /, q,) , (31) 
Fig. I2E: 

^ T^^-.-(p'; g^, / + k) S,{p - /)A„(/ + k) {-ieQ'^{2l + A;)^)A^(/) T'='".-(p; /, g,) , (32) 



Fig. EF: 



^2 ^ r^b,,c^^^,. ^^(^ _ i-teQ''^,)S,ip - l)AUl) T'=-'-(p; /, g,) , (33) 



with p' = p + k = pi + qi + k = pf + qf. 

For general momenta p, g (i.e. not necessarily on-shell), the quantity Aj^^'"* satisfies the 
relation 

k^Af^' = e| {Q' - Q^)A'^'''\q,p) - gM^^'"*(g,p + A;) + QM^^'"*(g - A;,p)} , (34) 
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Figure 2: Bubble chain diagrams for the process 705 — > (pB. Sohd, dashed and wavy hues 
correspond to baryons, mesons and photons, respectively. The square denotes the vertex from 
the leading order Lagrangian, the filled circle represents the bubble chain derived from the BSE. 

where A is the amplitude deduced from the Weinberg- Tomozawa term 

A'^'^\<1,P) = -^(^ + ^)<A^n[A^^A^A1> . (35) 

Making extensively use of Eq. ()34|) and the BSE it is straightforward to show that the 
contributions of Eqs. multiplied by fc^ yield in total 

el-Q-T'^'-^p'; qf, q,) +Q'T'^'-\p- q^, q^) -Q^T'i--\p'- g^, q^ + k) +Q^T'i-'^\p- qf-k, q,)] . (36) 

This compensates exactly the contributions from the remaining four diagrams where the photon 
couples to the external on-shell legs of the bubble chain. We have thus confirmed that gauge 
invariance is achieved if all possible diagrams of a photon coupling to a bubble chain are taken 
into account. In particular, it is not sufficient to consider only the coupling of the photon to 
external legs, since this will lead to a gauge dependent amplitude. On the other hand, within 
the field theoretical framework applied here it is also not sufficient to take into account only 
the coupling of the photon to the interaction kernel. Note that in the proof given here we do 
not assume the so-called on-shell approximation for the interaction kernel. 

It is also important to stress that an explicit evaluation of the BSE was not necessary 
and thus we do not need to specify the regularization scheme to render the loop integral in 
the BSE finite. Any regularization procedure which satisfies the Ward identities both for the 
propagators, Eq.(23), and the interaction kernel, Eq.(34), will maintain gauge invariance in the 
BSE as outlined in the proof. 
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Finally, the unitarity constraint for the full amplitude can be proven in analogy to the 
treatment in Sec. El The detailed derivation of unitarity is deferred to Sec. El Thus, also in 
order to obtain a unitarized amplitude (up to radiative corrections) one must take into account 
the coupling of the photon at all possible places in the bubble chain. 

4 Higher order interaction kernels 

In this section, we will consider more complicated structures for the interaction kernel as they 
arise at higher chiral orders in the effective Lagrangian. Since we restrict ourselves to cfP'BB 
contact interactions, the most general form of a term without the chiral invariant field strength 
tensor is given by 

Unt = B C,,...,,,,^,...,^,^^,...,^ [d^^ . . . [d^^^- . . . D^-0t j . . . d,.b) , (37) 

where we have suppressed flavor indices for brevity (we merely kept the symbol " ~ " as a 
reminder of the flavor structure indicating that the in- and outgoing baryons and mesons can 
be different) and introduced the notation = + HA^ with e the charge of the particle 
is acting on. As the contact interactions originate from a gauge invariant Lagrangian, charge 
conservation is guaranteed at each vertex. The introduction of explicit flavor indices does not 
change any of the following conclusions. Note that the constant C in Eq. (jSTj) niay also contain 
elements of the Clifford algebra. From this Lagrangian one derives both a (f){q)B{p) —>■ (f){q)B{p) 
contact interaction A{p,q,q) and a ■y{k)(f){q)B{p) — > (j){q)B{p) vertex e^A^{p,q,q,k), where 
is the polarization vector of the photon. Due to the form of the contact term (jHTj) . which 
can always be obtained by partial integration, the vertices do not depend explicitly on the 
momentum p. In Appendix El it is shown that they satisfy the relation 

/c^A^(pi, gj, g/, k) = e^A{pi, qi, qj - k) - e^A{pi, qi + k, qj) 

+ eBA{pi, qi, qf) - e^A{pi + fc, g^, g/) (38) 

for general momenta of the particles. This equation is the analog of Eq. ()34j) for the Weinberg- 
Tomozawa term. It follows then immediately by applying the same arguments as in the previous 
section that the coupling of the photon to the two-particle state of the BSE yields a gauge 
invariant amplitude also in the presence of more complicated contact interactions of the type 
Eq. dSZj). 

Note that throughout we have not considered the dimension two magnetic moment coupling 
~ cr'^'^f^^ and higher order operators involving the chiral covariant field strength tensor 
Such terms are of course present in the effective Lagrangian and must be considered at the 
appropriate order in the chiral expansion of the interaction kernel. However, these are of the 
form dpV^{0^^ — O^^) with some operator O^^ and the pertinent vertex in momentum space 
vanishes upon contraction with the photon momentum k^. 

As already mentioned in the previous section, the proof of gauge invariance does not de- 
pend on the specific choice of the meson and baryon propagators, but is rather valid for all 
propagators satisfying the Ward-Takahashi identities with the corresponding and '-yBB 
three-point functions. For example, one can define the BSE by employing propagators with the 
physical masses for the intermediate states and derive the interaction kernel from the effective 
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Lagrangian to a given chiral order. This automatically produces the correct physical thresholds 
of the unitarity cuts and we have followed this path in the present investigation. 

Alternatively, one may prefer to deduce the propagators from the effective Lagrangian as 
well. To leading chiral order this implies a common baryon octet mass shifting the threshold 
of the unitarity cuts to unphysical values. At higher chiral orders the inclusion of self-energy 
diagrams for the meson and baryon propagators will cure the situation by restoring the physical 
thresholds. The self-energy diagrams will modify the simple form of the propagators given in 
Eq. (OUj) : e.g., the meson propagators will acquire the form 

with M being the physical meson mass, Z the appropriate wave function renormalization 
constant, and the renormalized self-energy. 

In order to prove gauge invariance in the latter approach, one must also take into account 
the coupling of the photon to the self-energy corrections of the propagators. Since the effective 
Lagrangian is gauge invariant, the corresponding propagators and three-point functions satisfy 
the pertinent Ward-Takahashi identities and the proof is equivalent to the one given in the 
previous section. 

We also would like to compare the present investigation with the work of [11, 12]. Although 
similar in spirit, the authors study therein an integral equation for the two-body Green's func- 
tion, whereas we prefer to work with an equation for the scattering amplitude. "Gauging" the 
integral equation for the Green's function as outlined in [11,12], i.e. adding a vector index 
/X to all the terms of the equation such that a linear equation in /i-labeled quantities results, 
amounts to attaching an external photon everywhere including the external legs. In our frame- 
work, the straightforward application of the gauging method to the integral equation for the 
scattering amplitude fails, as in this case the procedure does not yield the contributions where 
the photon couples to the external legs. (Of course, one could correct this by adding the missing 
contributions by hand.) On the other hand, our approach is more convenient within the chiral 
effective framework, as it allows a direct comparison with the scattering amplitude derived in 
the perturbative scheme of ChPT. Moreover, the authors of [11,12] restrict themselves merely 
to one two-body channel, whereas in the present investigation this is generalized to several 
coupled channels. In contrast to [11,12] we explicitly specify the interaction kernel by deriving 
the vertices from the chiral effective Lagrangian and utilizing them as interaction kernels in the 
BSE. 



5 Unitarity 

After having constructed a gauge invariant amplitude for S07 Bcf) with Weinberg- Tomozawa 
or more general contact interaction kernels, we would like to investigate unitarity of the obtained 
amplitude. The calculation presented in this section generalizes the findings for the scalar field 
theory presented at the end of Sec. |2l as one must take care of the non- commutative nature of 
the matrix amplitudes due to the Clifford algebra and coupled channels. Moreover, we do not 
assume symmetry of the transition amplitudes under exchange of initial and final states. 
In operator form the statement of a unitary scattering matrix amounts to 

T - = -iT^T . (40) 
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For brevity we introduce a short-hand notation for the BSE, Eq. (j2EI), 

Tip) = A + J T{p)G{p)A = A + J AGip)Tip), (41) 

where p is the external momentum and G = iSA. The BSE for the meson-baryon scattering 
amphtude T is easily transformed into the unitarity relation 

T-f = j f{G - G)T (42) 

with O = 7o0^7o, as both T and G are elements of the Clifford algebra. Note that the adjoint 
also implies taking the transposed matrix in channel space. We will see that the quantity 
G — G is equal to setting the intermediate meson-baryon pairs on-shell in Eq. (jl^ . For invariant 
energies below the lowest three-particle threshold Eq. is thus equivalent to the unitarity 
constraint (jlOJ. 

If T is an analytic function, one can apply the residue theorem and rewrite the difference 
G — G as 

liSip - l)A{l) + S{p - l)A{l)) i{-27iif5+{P - M^)5+{{p - if - m') - / + m], (43) 

where I and p are the loop and external momentum, respectively, and 5+{k'^ — yU^) = 5{k'^ — 
li^)6{kQ). For a detailed derivation of this replacement, see e.g. [23]. The last equation can even 
be generalized to non-analytic T which does not have coinciding poles with G. In particular, the 
above replacement is valid if T describes the solution of the BSE with polynomial interaction 
kernels, as can be seen by insertion of its defining equation ()2(i|l into Eq. ()42j) . 

In order to prove unitarity for the transition 507 B(f), it is convenient to introduce the 
amplitude 

M^, = VLc + (44) 

with the disconnected piece 

and the transition amplitude calculated in Sections El and El The energies of the particles are 



given by Eg. = y Qi + and Ep- = \/ pi ^ m%, respectively. The piece V^^^^ is represented by 

the two disconnected diagrams in which the photon couples either to the baryon or meson, while 
the other particle does not interact at all. Although V^^^^ does not contribute to on-shell matrix 
elements and could in principle be omitted in the unitarity relation, its introduction generalizes 
unitarity beyond the physical region. From the amplitude M^^ one constructs the reversed 
amplitude M^^ for the process B{pf)(f){qf) B{pi)(f){qi)'-^{k). Neglecting radiative corrections 
and below the lowest three-particle (i.e. baryon two-meson) threshold unitarity implies 

- M>;^ = J f{p'){G{p') - G{p'))M^^ + J M>;^{G{p) - G{p))T{p) , (46) 

where we have replaced again the two-body phase space integration by the four dimensional 
integral over G — G. Note that in the physical region M^^ reduces to T^. By inserting the 
amplitudes and T from the BSE and making use of the unitarity statement for T, Eq. 
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as well as V^^^ = Vj^j^ and = from Eq. (j2Hl), one can indeed confirm the unitarity 
constraint for M^^ and thus for for on-shell matrix elements. 

We refrain from presenting the entire and tedious calculation here, but would like to com- 
ment on two points. First, the contribution of the disconnected graphs drops out on the l.h.s. 
of the unitarity statement ()46|) (due to the symmetry of these graphs under interchange of 
incoming and outgoing particles and V^/b ~ ^^/b)- r.h.s., they produce terms of the 

type 



Tip') iGip') - Gip')) 2E,X2n)'6^'\q, - l)Vi; 

= jfip') ii-2mf 5+iip' - I)' - m') 5+(/2 - M') [f'-t + m] 2E,X27r)'S^'Hq, - l)V^ 

= j dH fip')i-2ni) 5+iip' - /)2 - m") [/ - / + m] 5^'\q, - l)V^ 
= fip')i-2m) 5+iip' - q,)' - m')[f -^^ + m] 

= Tip') [Sip, + k)- Sip, + k)] VI (47) 

For on-shell matrix elements (and A; 7^ 0) this vanishes as e ^ in the propagators, but it 
happens that all terms of this type cancel each other on the r.h.s. of Eq. ()4fj|l even for general 
external momenta. Our second comment concerns some contributions from diagrams EE and 
F. On the l.h.s. of Eq. ()46|) one obtains, e.g., the combination 

I [ Tip') Sip' - I) Ail) Sip - I) Tip) + I [ Tip') Sip' - /) A(/) Sip - I) Tip) (48) 
Jl Jl 

which represents the discontinuity of the transition amplitude stemming from Fig. EF (the two 
integrals only differ in the sign of the "ie" terms in the propagators). According to the Cutkosky 
cutting rules and for momenta k"^ < 4m^ this discontinuity is given by 

Tip') iGip') - Gip')) Sip - /) Tip) + ^ Tip') Sip' - I) iGip) - Gip)) Tip). (49) 

In the kinematical region which is of relevance here the cut through the two propagators con- 
necting to the photon does not contribute to the discontinuity and can be safely neglected. We 
conclude by emphasizing that the unitarity constraint (j^6|) is only fulfilled if the photon couples 
to all possible places in the bubble chain. 



6 Conclusions 

In the present work, we have studied how gauge invariance is obtained for a photon coupling to 
a two-body state described by the solution of the Bethe-Salpeter equation. We have discussed 
the procedure both for a simple complex scalar field theory and for interaction kernels derived 
from chiral effective Lagrangians in the meson-baryon sector. In the latter case, we have first 
considered the Weinberg- Tomozawa term and afterwards the most general contact interaction 
consisting of two mesons and two baryons which can arise in the chiral effective framework. Our 
study underlines that it is not sufficient to take into account only the coupling of the photon 
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to the external legs as has been done in many calculations based on chiral unitary approaches, 
but one rather has to include all possible contributions of the photon coupling to the vertices 
and intermediate states. Neither is it sufficient to consider only the coupling of the photon to 
the interaction kernel. At the same time, coupling of the photon to the bubble chain at all 
possible places is necessary, in order to guarantee a unitary scattering matrix up to radiative 
corrections. 

For the interaction kernels discussed in the present work we have shown explicitly that 
gauge invariance is maintained in this manner. It is accomplished without assuming the on- 
shell approximation for the interaction kernel. Moreover, the explicit evaluation of the loop 
integral in the Bethe-Salpeter equation is not necessary and hence the proof does not depend 
on the chosen regularization scheme. But the regularization procedure is required to satisfy the 
Ward identities both for the propagators and the interaction kernels. This study will also be of 
importance for photo- and electroproduction processes of mesons on nucleons and for radiative 
decays of baryons and mesons which must be treated in a similar way, in order to achieve gauge 
invariant and unitarized amplitudes. Work along these lines is in progress. 
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A Higher chiral orders 

In this appendix we derive Eq. ()H8|) which follows from a contact interaction of the type 



A*l'"M;A'i + l'"MmMm + l'"A'n 



D^'K.. Df"^) ( D^'+i . . . D^'"0M ( D^'^+i . . . ) (A.l) 



with = + ieA^. The 'y{k)(f){qi)B{pi) (f){qf)B{pf) vertex is obtained from this La- 
grangian by extracting the part linear in Afj,. A partial derivative acting, e.g., on an incoming 
meson 0(g) yields a factor —iq in momentum space, while a partial derivative on (j){q)A{k) 
(both momenta incoming) leads to —i{q + k). Therefore, one obtains the vertex 



s=l 
m 

s=l+l 
n 

+ ^ J2 C, 



\(2«-2m+n) 



\(2«-2m+n) 



s-1 



=1 



lift" 



t=s+l 



.t=l+l 



t=l 
I 



s-1 



n (^/ - ^)"' 



.t=i+i 



+1 

rn 



t=m+l 



.t=s+l 



n p 

n 

n p 

t=m+l 



(-0 



(2«-2m+rt) 



s=m+l 



n 

.t=i 



.t=i+i 



s-1 

n 

J=m+1 



Up 

J=s+1 

(A.2) 
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Multiplying this equation with /c^ and making use of charge conservation at the vertex which 
follows from gauge invariance of the Lagrangian one arrives at 

- ^C',,..,„(-z)(^'-^"+")(-e^)gr ■ ■ ■ q^'iqj - kT^' ■ ■ ■ (Qf - kr-p^^' ■ ■ -pf" 

On the other hand, the vertex (l){qi)B{pi) — > (f){qf)B{pf) from the Lagrangian ()A.1|1 is given by 
the piece without the photon field A,_i and reads 

A{p,, g„ qj) = ■ ■ ■ g^'^^ ■ ■ ■ g^-pf-^^ ■ ■ -pf" . (A.4) 

This proves the Ward identity (jHHj) 

k^A^{pi, qi, qf, k) = e^A{pi, qi, qf - k) - e^A{pi, qi + k, qj) 

+ eBA{pi, qi, qf) - e^A{pi + fc, g^, qj) . (A. 5) 
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